<p> Draghilev Jacobian/rkfixed Solver </p>

44f7730c-b717-45ed-9552-3c8979546345 21 6 5 decimal

&[DATE] &[TIME] - &[FILENAME]

The Human eye contour project by Draghilev method

X0 tmin tmax N D n f rows 1 + : k 1 n range u k el x k el : for b 1 - (f u Jacobian 1 n 1 - n n submatrix * : матрица свободных членов ≡ A f u Jacobian 1 n 1 - 1 n 1 - submatrix : основная(квадратная) матрица ≡ n 2 > det 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : определителu дополнит.матриц ≡ det n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : det n el A det - : k 2 : k n 2 - ≤ k k 1 + : det k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : for 41 line det 1 el b : det n el A det : 21 line if dS 1 : i 1 n range D2 i el det i el dS / : for t x D1 D2 : Матрица  правых частей  системы ОДУ ≡ X0 tmin tmax N t x D1 rkfixed 91 line :

2.The parametrization "Human eye"

α 0.842 :

β 1.75 :

x y f1 1 x 2^y α - (y α + (* + (2^2 x * y α - (* (2^+ / β - :

Kоординаты начальных точек

R 0 y f1 0 ≡ y solve maple : 1.21033 1.21033 - 0.21671 i * 0.21671 i * - 41 mat

init1 0 R 3 el 21 mat : 0 0.217 i * 21 mat

init2 0 R 4 el 21 mat : 0 0.217 i * - 21 mat

f 1 el 1 x 1 el (2^x 2 el α - (x 2 el α + (* + (2^2 x 1 el * x 2 el α - (* (2^+ / β - :

Contours are defined by β ... solutions are hyper-reflexive.
'N' defines the fineness ... be generous.
Each contour closes on ≠ length of respective solutions
The best closed contours depend upon 't1max, t2max'. 
Adjust these two manually from the dual observer. tmin 0 : t1max 0.17977 : N 200 : t2max 0.0533 : 41 sys

B1 init1 tmin t1max N D :

B2 init2 tmin t2max N D :

dual observers wrt t1max, t2max B1 N row B2 N row 21 mat 0.178871 0.000097 - 0.433428 13 mat 0.053034 0.000203 - 0.433428 - 13 mat 21 mat

The parametrization Human eye by Draghilev method

B1 2 col B1 3 col augment B2 2 col B2 3 col augment 2 1 sys f1 B1 2 col B1 3 col augment B2 2 col B2 3 col augment 3 1 sys

ideally '0'

Both contours run clockwise

init1 0 0.217 i * 21 mat

init2 0 0.217 i * - 21 mat

k t x1 B1 k 1 + col t 1 + el eval :

k t x2 B2 k 1 + col t 1 + el eval :

τ 0 N 1 - range :